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Abstract. In the present paper a review and numerical comparison of 
a special class of multi-phase traffic theories based on microscopic, ki- 
netic and macroscopic traffic models is given. Macroscopic traffic equa- 
tions with multi- valued fundamental diagrams are derived from different 
microscopic and kinetic models. Numerical experiments show similari- 
ties and differences of the models, in particular, for the appearance and 
structure of stop and go waves for highway traffic in dense situations. 
For all models, but one, phase transitions can appear near bottlenecks 
depending on the local density and velocity of the flow. 



1. Introduction 

Traffic flow modeling has been considered on different levels of description, 
see [5] for a recent review: on the microscopic level the motion of each vehicle 
is described. Mathematical models are generally stated using a large system 
of ordinary differential equations for position and velocity of the vehicles 
based on Newtonian mechanics jTU HTJ El 1201 SO] . On the macroscopic level 
the state of the system is described by averaged quantities. Typically, den- 
sity and linear momentum are used to describe the flow. The corresponding 
mathematical models are based on systems of nonlinear partial differential 
equations derived from conservation laws with suitable closure relations. 
Starting from the pioneering work of Aw and Rascle [2] new macroscopic 
models for traffic flow have been derived and investigated intensively in the 
last decade, see for example HH HJ [23l EH [28] . These models avoid several 
inconsistencies of previous models, like wrong way traffic and missing bounds 
on the density. We note that these models can be derived from microscopic 
models in a variety of ways, see for example [TJ \52\. Finally, kinetic theory 
describes the state of the system by a probability distribution function of 
the position and velocity of the vehicles. Mathematical models generally use 
integro-differential or Fokker-Planck type equations. Kinetic equations for 
vehicular traffic can be found, for example, in |45 | I44 | l4"2j 138]. Procedures 
to derive macroscopic traffic equations including the Aw/Rascle model from 
underlying kinetic models have been performed in different ways by several 
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authors, see, for example, [27] and [39]. These procedures are developed in 
analogy to the transition from the kinetic theory of gases to continuum gas 
dynamics. 

Another basic problem of macroscopic traffic flow equations has been 
described by Kerner [S2J [S31 [33]. The observations there suggest a more 
complicated dependence of the homogeneous steady speed states on den- 
sity: these states are not given by a uniquely defined function u = U e (p) 
as in the above mentioned models, but cover a whole range in the density- 
flow diagram leading to a multi-valued fundamental diagram. The resulting 
dynamical system has a multi-phase behavior in the sense of Kerner: the 
flow changes between different stationary state which represent free and so 
called synchronized or jam behavior. In the context of the derivation of 
macroscopic models from microscopic ones the homogeneous steady state 
solutions can be interpreted as an emergent behavior of interactions at the 
microscopic scale and multiple solutions may be related to the heteroge- 
neous behavior of the driver-vehicle subsystem. A variety of microscopic 
and macroscopic multiphase models has been developed by several authors. 
In particular, there is a large number of works on microsopic models. We 
refer among many others to [361 [371 S3 SZ1 EH [53J EH EI] . Macroscopic 
models for traffic flow with phase transitions in the sense of Kerner can be 
found in EU [T3j [12] . However, these models do not describe phenomena 
like stop and go waves near bottlenecks. For microscopic and macroscopic 
multiphase models exhibiting stop and go waves and similar traffic instabil- 
ities we refer to [361 E3 E31 HOI E2J HE] . Kinetic models allowing for multiple 
stationary solutions and associated macroscopic models with multi-valued 
fundamental diagrams have been developed in (29] [32j [26j 03] . We refer to 
[5j for a recent review and to [17} HI [16] for further material on the above 
issues. 

The present paper contains a comparison and discussion of a class of 
macroscopic models with multi- valued fundamental diagrams. We consider 
models of the form 

(1) d tP + d x (pu) = 0, 

dt(pu) + d x (pu 2 ) - c(p)d x u = pR(u, r) 
with right hand side 

R(p,u) = £ [U(p,u)-u] 

and fundamental diagrams given by functions U = U(p, u) having at least 
two equilibrium solutions, i.e. solutions of the equation u = U(p, u) for fixed 
p out of a certain density domain, where multiphase traffic may appear. 

The paper starts with a review of this class of macroscopic multi-phase 
models. The models are derived from either microscopic or kinetic equa- 
tions. To guarantee a proper comparison of the considered models several 
changes to the original models are proposed. Moreover, the parameters of 
the different models are chosen such that the stable equilibrium solutions are 
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the same for all models. Then, the different models are numerically investi- 
gated for a bottleneck problem and the appearance of stable wave patterns 
is shown which can be interpreted as stop-and go waves at (on ramp) bot- 
tlenecks. This numerical comparison as well as the changes made for each 
of the models to make them comparable are new up to the knowledge of 
the authors. We note that stable periodic waves excited by small periodic 
perturbations have been studied in a series of papers also for equations with 
single valued right hand sides, see |24 | 1251 ITU] . Remarks on these waves can 
be found in section [5j 

The paper is arranged in the following way: In Section [2] the derivation of 
macroscopic equations from microscopic models is reviewed and applied to 
a multi- phase traffic model from |35| [36] . In Section [3] kinetic equations are 
investigated and used to derive multi-valued fundamental diagrams. The dif- 
ferent models are partially changed to make them comparable to each other. 
Section [4] contains a summary and comparison of the different approaches 
and the derived multi-valued fundamental diagrams. Finally, in Section [5] 
numerical results are given comparing the different density-velocity rela- 
tions. Moreover, an inhomogeneous traffic flow situation with a bottleneck 
is investigated, showing the appearance of traffic instabilities together with 
a qualitative comparison of the structure of these instabilities. 

2. Continuum multi-phase traffic model derived from 
microscopic equations 

2.1. From microscopic to macroscopic models. We review the classical 
procedure for so called 'General Motors' (GM) type car-following models, 
see [HI HI]. Denoting with Xi(t) , Vi(t) , i = 1, . . . , N the location and speed 
of the vehicles at time t G M + , and the distance between successive cars by 

we consider the microscopic equations 

X{ — Vi 

(V i+ l - Vi) 1 

Vi = C { _ R +-(U(pi)-Vi). 

The local "density around vehicle i" and its inverse (the local (normalized) 
"specific volume") are respectively defined by 

H Ik 

Pi = — and n = — = — , 
k Pi ti 

where H is the length of a car. 

Remark 2.1. Here, the density is normalized and therefore dimensionless, 
so that the maximal density is p m = l/r m = 1. 

The constant C > and the relaxation time T are given parameters. The 
function U = U(p),0 < p < p m — 1 is the so called fundamental diagram, 
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see [HI EH HE]. The simplest choice is given by U(p) = 1 — p. One obtains 
the microscopic model 

(2) ±i = v i , 

C (v i+1 -Vi) 1 



We have 



ii = v i+1 - Vi or n = jj (vi+i ~ Vi) 



The limit of number of cars going to infinity yields the Lagrangian form of 
the macroscopic equations, see p]. We obtain the equivalent of the p-system 
in gas dynamics (isentropic Euler equations in Lagrangian form), compare 
1221 , 



(3) Ott — dxu = , 

d T u - c{p)d x u = -[U (p) - u] , 

where r = - is the specific volume, i.e. the (local) dimensionless fraction 
of space occupied by the cars, p the (normalized) density is the limit of pi 
defined above, as the number of cars tends to infinity, u is the macroscopic 
velocity of the cars. Moreover, 

(4) c(p) =c(--l 

and the function U(p) is defined in the microscopic model above. We change 
the Lagrangian "mass" coordinates (X, T) into Eulerian coordinates (x, t) 
with 

d x x = P , d t x = - P v, t = t 

or 

dxx = p~ l = r, 8tx = v. 

Thus, X = f x p(y, t)dy describes the total space occupied by cars up to 
point x. The macroscopic system in Eulerian coordinates is then 

(5) d t p + d x {pu) = 0, 

d t (pu) + d x (pu 2 ) - c(p)d x u = ^ [U(p) - u] . 

For the well-posedness of the above problem we refer to |S]. 

Remark 2.2. The same approach works for right hand sides with multi- 
valued equilibrium distributions 

(6) | [U(p,u)-u]. 

Examples are Switching curve (SC) models as in [25] with 
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U(p,u) 



\Ui(p), p < pj oru> S(p),pf < p < pj, 
[U 2 (p), u < S(p),p f < p < pjor p > pj. 



Here the switching curve S = S(p) is given. We consider the density in 
the (synchronized flow) region between a lower bound of free flow pf and an 
upper bound of jam traffic pj. Then, there exists multiple stationary states 
U\,U2 whose regions of influence are separated by the switch-curve S. 

A similar model is the Speed- adaptation (SA) type models of Kerner et.al 
36 with 



U(p,u) 



\Ui(p), U > U sync , p < Pj , 
\U 2 (p), U <U sync ,p> pf, 



where the parameter U syn is the averaged speed, which separates the domains 
of influence of the two stationary states in the 2-D region of synchronized 
flow in the flow-density plane. 

2.2. A microscopic ATD-type model. In this section we sketch a sim- 
plified version of the microscopic Acceleration time delay (ATD) model of 
Kerner et. al [36]. We consider a microscopic model with the variables 
space, velocity and acceleration: 



(7) 



X j 
Vi 



a i 



with 



Avi 



F{ai,Vi, -jf, T i) 



F(ai,Vi, 
:a{ ree -a i )/T l 



V-i+l - Vi X i+ l - Xi 



H 



H 



syn 



jam 



del i 
Ui)/T d el, 

ai)/Tdeh 



n > G(vi) , T% > Tjam, 

, T{ > Tj am , 



Here, Avi = Vi + \ — Vi and a{ ree , a^ n , a{ am denote the desired accelerations 
in the free, synchronized and jam region respectively and T<j e z denotes the 
time delay of the acceleration of the vehicle. The function G separates the 
free from the synchronized acceleration behaviour and will be fixed later at 
the end of Section [4} For a proper comparison with the above models we 
change the definitions in|36| of the different accelerations slightly and define 
the terms as follows 



free 



jam 



1 



Vi) + 



Avi 



mm(U(pi) 



c(Pi, 
H 

v u 0) + ^Ar, 



H 



1 c(Pi) . 

- fVl + ^lAvi, 
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where U{p) and c{p) are given as before. This means acceleration depends 
on the speed difference to the predecessor and a term relaxing to a desired 
acceleration. 



The hydrodynamic multi-phase model. To obtain the hydrodynamic version 
of the microscopic model in the last section we follow the procedure in 
Section 2.1 In Lagrangian coordinates we obtain directly 

(8) d T r - d x u = , 

&tu = a 

dxa = F(a,u,dxu,r). 
This leads to the following equations in Eulerian coordinates 

(9) d tP + d x (pu) = o, 

dt(pu) + d x (pu 2 ) = pa 

d t (pa) + d x (pua) = pF(a,u,Td x u,r). 

A reduced model. Assuming that the delay times for acceleration are small 
we can reduce the above ATD-type model to 



(10) 



where 



d t p + d x (pu) = 0, 
d t (pu) + d x (pu 2 ) - c(p)d x u = pR(u, r), 



(1 (U(p)-u 



R(p, u) 



y V ^ VF7 — u,j, T > G(u),T > Tj, 

y min(C/(p) — u, 0), r < G(u), r > Tj, 



rp LI, 



T < T,' 



This is equivalent to 



R(p, u) 



Mp,u) 



u 



with 



U{p,u) 



U(p), t > G(u), p < pj or t < G(u),u > U(p), p < pj, 
u, t < G(u),u < U(p),p < Pj, 
0, p> pj. 



For comparison with the other multi-valued fundamental diagrams we 
rewrite the relaxation term using K(u) = 1/G{u): 



U(p,u) 



U{p), p < K(u),p < pj or p> K(u),u> U(p),p f < p < pj, 
u, p> K(u),u <U(p),p f < p < pj, 
0, p > pj. 



Continuum multi-phase traffic theories 



7 



3. Multi-phase hydro-dynamic equations derived from kinetic 

equations 

3.1. Kinetic equations and correlations. The basic quantity in a kinetic 
approach is the single car distribution f(x, v) describing the density of cars 
at x with velocity v. The total density p on the highway is 

P( x ) = / f(x,v)dv, 
Jo 

where w denotes the maximal velocity. Let F(x, v) denote the probability 
distribution in v of cars at x, i.e. f(x,v) = p(x)F(x,v). The mean velocity 
is 

/•w 

u(x) = I vF(x,v)dv. 
Jo 

An important role is played by the distribution f( 2 \x,v,h,v+) of pairs 
of cars being at the spatial point x with velocity v and leading cars at x + h 
with velocity v+. This distribution function has to be approximated by the 
one-vehicle distribution function f(x,v). We use the chaos assumption 

/ (2) (x, v, h, v + ) = q(h) f(x, v) F(x + h, v + ), 

compare |42|. For a vehicle with velocity v the function q(h, v; p, u) denotes 
the distribution of leading vehicles with distance h under the assumption 
that the velocities of the vehicles are distributed according to the distribution 
function /. 

Thresholds for braking (Hb) and acceleration {Ha) are introduced. From 
a microscopic point of view drivers will brake, once the distance between the 
driver and its leading car is becoming smaller than a threshold Hb and will 
accelerate, once this distance is becoming larger than Ha- Otherwise the 
cars will not change the velocities. Velocities are changed instantaneously 
once acceleration or braking lines are reached. Models including acceleration 
of the cars can be developed as well, see [32] for an example. 

The distribution of leading vehicles q(h) is prescribed a priori. The main 
properties, which q(h) has to fulfill are positivity, 



POO 

/ q{h)dh = 1, 
Jo 



and 

1 



(11) / hq(h)dh 

P 



Equation (11) means that the average headway of the cars is 1/p. The 
leading vehicles are assumed to be distributed in an uncorrelated way with 
a minimal distance Hb from the car under consideration, see 

q(h) = pe-^ h - H ^ X[HBtOQ) (h). 
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The reduced density p has to be defined in such a way, that ( 11 ) is fulfilled. 
One obtains 

We note that 

q A = q(H A ) = pe-K H *- H ») 

and 

q B = q(H B ) = p ■ 

The probability P ov = P ov (p,u) for overtaking or lane changing and the 
corresponding probability Pg = 1 — P ov for braking are determined from 
phenomenological considerations: at constant density, free flow of cars, i.e. 
larger velocities will be related to larger probabilites of overtaking or smaller 
probabilites of braking. So called synchronized traffic is associated to smaller 
velocities and thus larger probabilites of braking. That means the probabli- 
tiy of braking can be considered as - for fixed density - a decaying function 
of velocity u. Similar arguments can be found for example in [35J. 

Remark 3.1. In the following we present a kinetic model. Note that the 
results like multi-valued fundamental diagrams and stop and go behaviour of 
the derived macroscopic equations do not depend on the exact choice of the 
microscopic interactions we have chosen here. The model discussed in the 
next section is only chosen due to the fact that explicit stationary solutions 
are available. We could as well have chosen models like in [39] or Fokker- 
Planck type models like in |32j . 



3.2. The evolution equation. To write the kinetic evolution equations in 
a simple form we use 

k = k(p, u) 



q A + PsqB 
and 

7 = l(p,u) = = qA + PBqB- 

We consider a relaxation frequency 

v = v(k) 

and define 

1 

The kinetic model is then given by the following evolution equation for 
the distribution function /: 
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(13) d t f + vd x f = C + (f) 

k(G+-L+)(f) + (l-k)(G+-L+)(f) 

+u(G s -L s )(f) 
with the loss and gain terms for braking interactions 

G bU) = \v-v + \a B (v;v,v + )f(x,v)F(x + H B ,v + )dvdv^ 

J Ji)>VJ l . 

L+(f) = [ \v-v + \f(x,v)F(x + H B ,v + )dv + . 
J v+<v 

The loss and gain terms for acceleration is defined as 

G \W) = \v-v + \a A (v]v,v + )f(x,v)F(x + H A ,v + )dvdvj, 

J 

L+(f) = f \v-v + \f(x,v)F(x + H A ,v+)dv+. 
Finally terms describing the random behavior of drivers are 



G s{f) = / a s (v,v)f(x,v)dv 
Jo 

Ls(f) = f(v). 

as and a a denote the distribution of the new velocities v after an in- 
teraction. Reaching the braking line the vehicle brakes, such that the new 
velocity v is distributed with a distribution function as depending on the 
old velocities v,v+. For acceleration the new velocity is distributed accord- 
ing to a a- The relaxation term is introduced to include a random behaviour 
of the drivers. 

Remark 3.2. For further details on this Boltzmann/Enskog approach to 
traffic flow modelling see [38\ . 

Example. For the probability distributions a a, o~b we choose the following 
simple expressions: 

(14) <T B (V,V,V+) = ; \ X{v+,v](v) 

and 

(15) a A (v,v,v+) = - 1 x X[v,v+](v)- 

This means we have an equidistribution of the new velocities between the 
velocity of the car and the velocity of its leading car. Finally, 

(16) <rs(v,v) = — . 

w 
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3.3. Stationary Distributions and multi-valued Fundamental Dia- 
grams. In this section we investigate the stationary homogeneous equations 
and determine the multi- valued fundamental diagrams. We consider the lo- 
cal interaction operator: 

C(f) = 7 [k(G B - L B )(f) + (1 - k)(G A - L A )(f) + u{G s - L s )(f)] 

with / = pF. The gain and loss terms Gb,Lb, etc. are defined in the same 
way as G^,L^, etc., except that x + Hx,X = A, B is substituted by x, 
wherever it appears. The homogeneous stationary equation is 

C(f) = 0. 

We assume that for fixed p and k there is a unique solution 

f = r = P F £ (k,v) 

of this equation. This is true for the example stated above. 
Thus, for fixed k the mean value of F e is then 

rW 

u e (k) := / vF e (k,v)dv. 
Jo 

The function u e is uniquely determined due to the above assumption as a 
function of k. However, this does not yield immediately the fundamental 
diagram, i.e. an equilibrium relation between flux and density. 

Instead, the fundamental diagram is determined from the following con- 
siderations: let u be the (possibly multi-valued) solution of the equation 

(17) u = u e (k(p,u)) 

for fixed p. If there is a unique solution we obtain a well defined relation 
for equilibrium velocity and density and the usual fundamental diagram. 
However, in general this equation will have a multitude of different solutions 
u, even infinitely many. Plotting a dependence of this solution on the density 
one obtains in the general case a two-dimensional region in the density- 
velocity plane, where the solutions are located. The fundamental diagram 
is then a multi-valued function. 

Remark 3.3. For the example above the homogeneous solution can be solved 
explicitly and the corresponding multi-valued solutions of equation H| can 
be evaluated numerically. Explicit expressions for F e {k) and u e (k) can be 
found in |26j. A plot of u e (k) is shown in Figure^ 

Remark 3.4. In contrast to the other models described above the kinetic 
approach gives an explanation for a multi-valued fundamental diagram using 
in particular the braking probability Pb as a basic quantity. 

3.4. Derivation of Macroscopic Models. In this section macroscopic 
equations for density and mean velocity are derived. Different procedures 
are described, for example in [391 EH- 
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Balance Equations. Multiplying the inhomogeneous kinetic equation ( 13 ) 
with 1 and v and integrating it with respect to v one obtains the following 
set of balance equations: 

(is) d tP + d x {pu) = o 

dt(pu) + d x (P + pu 2 ) + E = S 
with the 'traffic pressure' 

(19) P = / (v- uffdv, 



the Enskog flux term 

/•w 

(20) E = / v[C(f)(x,v,t)-C+(f)(x,v,t)]dv, 

Jo 

and the source term 

j"U] 

(21) S = / vC(f)(x,v,t)dv. 



For the present discussion we are, in particular, interested in the source 
term S. 

Closure. We concentrate on the relaxation term and cite the results for the 
other terms , compare [39] ■ The traffic pressure P is negligible and approx- 
imated by zero, see [391 . Moreover, the Enskog term E is approximated by 



linearizing expression (20) for E in H. We obtain [39] 

E ~ -c k in{p)d x u 

with Ckin{p) given by the details of the collision operator. In the following 
we will neglect the special form of Ckin and choose Ckin = c(p) for comparison 
as in the other models described above. 

Finally, the source term S has to be approximated. We use a relaxation 
time approximation 

c(f) ~ i (nk( P ,u), v)-f( V )). 



This yields 



S ~ S e (p,u) = p±( u t(k(p,u))-u). 



Thus, from the kinetic approach one obtains macroscopic equations of the 
form 

(22) d t p + d x {pu) = 

dt(pu) + d x (pu 2 ) - c(p)d x u = S e (p, u) 

with 

S e (p,u) = p^(u e (k(p,u))-u), 
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where k = k(p, u) is defined as 

PbQb 



k 



QA + PbQb 1 + exp(-p(H A -H B ) 



Choosing Ha = Hb this simplifies to 

k 1 



Remark 3.5. One obtains a multi-valued variant of the Aw-Rascle equa- 
tions with a multi-valued relaxation term on the right hand side. 

4. Comparison of multi-phase hydrodynamic models 
We consider models of the form 

(23) d tP + d x (pu) = 0, 
d t (pu) + d x (pu 2 ) - c(p)d x u = pR(u, r) 

with 

(24) c(p) = c(--l 



and 



P 



R(p,u)= p - [U(p,u)-u] 



and fundamental diagrams given by functions U of the following form: 
ATD-type models. 

U(p), p < K(u),p < pj or p> K(u),u> U(p),p f < p < pj, 
U(p, u) = I u, p > K(u),u < U(p), pj < p < pj, 
0, P>Pj- 

SA-type models. 

\Ui(p), U > U S y n c,p < Pj, 

\u 2 (p) U <U sync ,p> p f , 
See Kerner et.al. [36] for details. 



U(p,u) 



Switching curve models. 

\Ui(p), p < p f or u> S(p),p f < p < pj, 
\u 2 (p), u < S(p),p f < p < pj or p> pj. 

where S(p) is a switching curve. For an investigation of these models, see 



U(p,u) 
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Kinetic models: 

U(p,u) = u e {k(p,u)) 

with 

l 

k 



In Figure [3] we plot the equilibrium solutions of u = U(p, u) together 
with the values U(p,u) — u denoting the length of the arrows. Moreover, 
U(p, u) — u is plotted for fixed p versus u in Figures [7| 

For a proper comparison of the above models the parameters are chosen 
as follows: 



U(p) = U 1 (p) 
K- 1 {p) = U 2 {p),p j > p> p f 
K(u) = Pf,u> U2{pf) 

U 2 {p) < U sync <Ui{p),Pf < p < Pj. 

Moreover, the functions v and Pb in the kinetic model are chosen such 
that the stable kinetic equilibrium solutions of u = u e (k(p,u)) are given by 
Ui and U2, the instable solution by S(p). This leads to the fundamental 
diagrams shown in Figures [3] to [6} 



5. Numerical Investigations 

In this section we compare the different approaches numerically. First, 
the kinetic model is investigated and the associated fundamental diagram 
is determined. Second the other three models stated in the last section are 
compared to the kinetic approach and third all three approaches are used 
in a inhomogeneous traffic simulation with a bottleneck. 



5.1. The stationary, homogeneous kinetic equation. We consider the 
kinetic equation and resulting fundamental diagrams. For the numerical 
simulations we normalize and use w = 1. 

Moreover, we choose v as in Figure [T] A reasonable function v should 
be zero for maximal density (k = 1). In this case there is no more random 
behaviour of the drivers, all drivers have velocity 0. For the case k = we 
have chosen v as a finite quantity. If these two features are fulfilled, the 
qualitative behaviour of the model does not depend on the exact form of v. 
The braking probability Pb is plotted as well in Figure [TJ 

Using Pb and v described above we compute for fixed k the unique sta- 
tionary solution of the homogeneous kinetic equation and the function u e (k) 
following Section 3^ The dependence of u e on k is plotted in Figure [2} 



11 



Raul Borsche,Mark Kimathi and Axel Klar 




Figure 2. Function u e (k) 

5.2. The multi-valued fundamental diagrams. In this subsection we 
plot the multi-valued fundamental diagrams for the four cases discussed in 
section |3j The functions U\(p) and ^(p) for ATD and switching curve 
models were chosen as 

Ux(p) = U Q Uvh(^- (--0.05 
\T U \p 

VM =[ /; t anh(^(i-l.l)) 

with Uq = 0.85, C v = 0.45, U£ = 0.5, T = 2.9. Moreover, T = 5, C = 
0.3, U sync = 0.28 and S = S(p) is given by a linear function connecting 
Ui(pf) and U2(pj)- The solutions of the nonlinear equation u = U(p,u) 
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Figure 3. U(p,u) — u for the kinetic model. On the Right: 
a zoom of the multi-valued region. 




Figure 4. U(p,u) — u for the switching curve model. On 
the Right: a zoom of the multi- valued region. 

are plotted together with a plot of the quantity U(p, u) — u as arrows with 
direction ±u. Figure [3j [4j [5] and [6] show the multi-valued fundamental 
diagram (speed-density relation) for the different models. In each figure, a 
zoom of the multi-valued region is shown. 

In all cases the values for pf and pj are chosen as 0.3 and 0.5 respectively 
such that for p < pf we have only one steady solution. For pf < p < pj 
three (infinitely many for ATD) solutions exist. And for the region p > pj 
again only one solution exists. 

Moreover, Figure [7] shows the values of U(p, u) — u for a fixed value p = p 
with pf < p < pj. 

5.3. Numerical solution of multi-phase macroscopic equations. Fi- 
nally, the macroscopic equations are investigated for a bottleneck situation. 
For the computations we choose a Godunov method, see |49j . We use a mesh 



16 



Raul Borsche,Mark Kimathi and Axel Klar 




Figure 5. U(p,u) — u for the SA model. On the Right: a 
zoom of the multi-valued region. 




Figure 6. U(p, u) — u for the ATD model. On the Right: 
a zoom of the multi-valued region. 



size Ax = 0.15, the Courant number X c fi = 0.99 and a computation time 



Tend = 400 units. Figures pi 10 and 11 show the velocity and density in 



space and time for a three lane highway with a reduction of lanes from 3 to 2 
at X = for the four different models. In the simulation, the lane reduction 
is achieved by multiplying the density p in the term on the right hand side of 
the equations by a factor | for X > units and using a linear interpolation 
between the 2 regions. Apart from the ATD-type model, one clearly observes 
large changes in velocity and density in the backwards propagating traffic 
jams which might be interpreted as stop and go behaviour. Figure 12 show 
the flow-density relation at various locations of the considered highway, i.e. 
upstream of the bottleneck (X = —20), within the bottleneck (X = 0),and 
downstream of the bottleneck (X = 5). The flow rate drops from the initial 
value used in the simulation to settle at the maximum values for the con- 
sidered highway's downstream location, X = 5. For the ATD model in its 
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Figure 7. U(p,u) 
and ATD models. 



u for the kinetic, switching curve, SA 
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Figure 8. Spatio-temporal congested traffic pattern - veloc- 
ity (Left) and density (Right) for the kinetic model. 



present form we obtain a rather different behavior due to the zero forcing 
inside the multi-valued region. 

Remark 5.1. The models differ in the frequency and uniformity of the 
waves generated at the bottleneck. The fact that the ATD model does not 
generate stable waves in this situation does not mean that the model is in 
general incapable to describe traffic situation with such patterns. For exam- 
ple, in [36] situations are described where these waves appear. The models 
derived from the kinetic equations can be viewed as intermediate models be- 
tween switching curve and ATD model, compare Figure^ 

Remark 5.2. A similar investigation could be performed for the model in 
[12] with suitable right hand side. 
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Figure 10. Spatio-temporal congested traffic pattern - ve- 
locity (Left) and density (Right) for the SA model. 
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Figure 11. Spatio-temporal congesting traffic pattern - ve- 
locity(Left) and density (Right) for the ATD model. 
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Figure 12. pu — p relation for the kinetic, switching curve, 
SA and ATD models upstream (x = —20), within (x = 0) 
and downstream (x = 5) of the bottleneck. 
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Remark 5.3. Using the coefficient c{p) = Cp as done for example in |14j . 
we obtain similar simulation results as above, if the parameters are suitably 
choosen. 

Remark 5.4. The stable waves excited by small periodic perturbations as 
discussed in |25} 124"] which may also be obtained from equations with single 
valued right hand sides are usually not persistent anymore for bottleneck 
situations. These waves are damped out as the high density region travels 
backward from the bottleneck. 

Summary. Multi-valued fundamental diagrams are obtained using different 
approaches: a derivation from microscopic equations given in [34, 36J, from 
kinetic models as in [26] and a phenomenological macroscopic model from 
|25j . These approaches are compared with each other from the point of 
view of their multi- valued fundamental diagrams and for an inhomogeneous 
bottleneck simulations without any external excitation. Apart from the 
ATD-type model, all the other models are able to show stop and go patterns 
for the described situation with a bottleneck without external excitation of 
waves by the ingoing flow from an on-ramp. 
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